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Abstract. The Hartree energy shift is calculated for a unitary Fermi gas. By 
including the momentum dependence of the scattering amplitude explicitly, the 
Hartree energy shift remains finite even at unitarity. Extending the theory also 
for spin-imbalanced systems allows calculation of polaron properties. The results 
arc in good agreement with more involved theories and experiments. 
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1. Introduction 

Strongly interacting atomic gases have proved to be an excellent platform for studying 
many-body quantum phenomena [JJ. The purity and the diluteness of the gases 
allow the description of the system using well understood microscopic models [H [31 0] 
and a host of theoretical tools have been applied for solving stationary and dynamic 
properties of the system [6l [3 [8] . 

Despite the multitude of more exact and involved theoretical tools, mean-field 
theories form an important class of theories due to their inherent simplicity and the 
applicability to three-dimensional systems [5]. Surprisingly, these mean-field theories 
have proved to be able to describe qualitatively even strongly interacting systems |10j 
although finite temperatures and polarized gases have proved challenging. Indeed, 
critical temperatures and critical polarizations for the superfluid-normal transitions 
are exaggerated in a mean-field theory. These problems are caused by the difficulty 
of describing the non-superfluid state of the gas properly, the state which is often 
referred to as the 'pseudogap phase' [5]. Indeed, the non-superfluid state in the BCS 
mean-field theory is a non-interacting finite temperature Fermi sphere in which all 
interaction effects are absent; a strong assumption for a strongly interacting system. 

The natural question is whether the simple BCS mean-field theory can be 
improved to include interaction effects in the normal state at least approximately. 
Obviously, there are a lot of extensions of the simple BCS theory [TT1 IT21 [T31 [HI 1151 [TBI 
OH [T71 HH1 [THl HOI HI] that incorporate fluctuations in the BCS model and can describe 
the non-superfluid state properly. Most of the methods are however numerically very 
heavy and the simplicity of the BCS mean-field theory has been lost. The approach in 
this work is quite different. Instead of considering fluctuations around the BCS theory, 
I neglect all pairing effects altogether and concentrate on one key approximation done 
in the BCS theory when in the strongly interacting regime, namely the omission of the 
Hartree energy shift. The idea is not new as pairing correlations are often neglected 
when considering systems at high temperatures but here I will relax approximations 
often done to the Hartree energy shift [22] • The result is a self-consistent theory which 
is simple but sufficiently powerful to describe a wide range of phenomena. 

The structure of this work is the following: in Section [2] I will first provide a 
definition for the Hartree energy shift for a strongly interacting gas. In Section [3] 
the main points of the standard scattering calculation for deriving the many-body 
scattering T-matrix are reviewed. The details of the derivation can be found 
in |Appcndix A| In Section [4] the Hartree self-energy will be calculated for a balanced 
two-component Fermi gas at unitarity. In Section [5] the theory is generalized for spin- 
imbalanced systems and the Hartree self-energy and the effective mass are calculated 
for a polaron. In Section [6] the theory will be further extended for calculating the 
occupation numbers at the high momentum limit and the Tan's contact coefficient 
will be solved both for a balanced system and for a polaron. The key approximations 
are reviewed and the implications of the theory are discussed in Section [7] Finally, 
the key findings are summarized in Section [8] 
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2. The Hartree shift 

The effect of atom-atom interactions enter the single-particle Green's function G(P) 
through the self-energy £(P) according to the Dyson equation 

G(P) = G (P) + G (P)E(P)G(P) - ^(p)-!,^) . (!) 

where the non-interacting Green's function G (P) = G {p,Pa) = l/(Po — e p)- The 
standard approximation for the self-energy used in the context of dilute atomic gases 
is the Hartree-Fock self-energy written using a scattering T-matrix 

E(P) = "* / t|ti G ( K ) r (K, P; P (2) 

This describes the self-energy of a particle with four-momentum P, interacting with 
a particle of four-momentum K drawn from the bath. The interaction is described by 
the many-body scattering T-matrix T(K, P; P, K) which promptly returns the particles 
back to the initial state, thus returning the particle with the four-momentum K to the 
bath and the initial particle to the state P. The self-energy utilises only the diagonal 
part of r, but the T-matrix will eventually be needed for a more general case as well. 
The above form assumes that the two interacting particles are distinguishable and 
thus it suffices for describing interspecies interactions in two-component Fermi gases. 
For indistinguishable particles one should add also the exchange interaction channel. 

Replacing the many-body scattering T-matrix by a constant V m _r a, where 
a is the s-wave scattering length, yields the well known result for the Hartree energy 
shift in the weakly interacting limit 

^Hartree — an i 

m 

where n is the atom density. However, this formula does not work in the strongly 
interacting, or unitary, limit where a — > oo. Indeed, at unitarity the only relevant 
length scale in a dilute Fermi gas is the Fermi momentum fcp and the Hartree energy 
shift must become independent of a [23|[TD]. Simply replacing a by 1/fcp would yield 
as the Hartree energy shift 4/3irEp 0.42 Pp, which is already in a rather good 
agreement with the values obtained for the Hartree-like energy shift in a Monte Carlo 
study [24] and a self-consistent Green's function theory [FT] . 

However, such a naive argument raises more questions than answers, and the 
purpose of this work is to make the argument rigorous and especially clarify the 
physics behind it. However, the main problem with discussing the Hartree shift in the 
strongly interacting regime or at unitarity is that there is no good definition for it. 
While in the weakly interacting limit the Hartree shift can be considered as an energy 
shift linear in density n, at unitarity it is proportional to n 2 / 3 . Another important 
issue to notice is that the Hartree shift is not constant, i.e. equal for all atoms, even 
in the weakly interacting limit. This follows from the breaking of the translational 
symmetry due to the presence of the Fermi sphere [35] coupled with the momentum 
dependence of the two-particle scattering amplitude, which essentially says that the 
interaction strength between two atoms becomes weaker when the relative momentum 
increases. For weakly interacting atoms this momentum dependence is not relevant 
as long as |fcpa| <C 1, but in strongly interacting systems one can hardly expect the 
Hartree shift to be a constant. 
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Since the Hartree shift in the strongly interacting systems has so few properties 
in common with the weakly interacting limit, one can ask whether it even makes sense 
to talk about Hartree shifts at unitarity. Indeed, many advanced methods designed 
for describing strongly interacting systems do not separate the Hartree contribution 
at all [TJJ [T5] or identify it only as a part of the total grand potential [35] , which is a 
global quantity and not a quantity directly relevant for single-particle excitations (and 
hence for the single-particle Green's function G{K)). However, in simpler mean- field 
theories, such as the BCS theory applied to the strongly interacting regime [371 [35] . 
the Hartree shift is explicitly neglected and consequently also easily identified [101 . 
The approach of this work is to omit the terms traditionally included in the BCS 
theory and instead consider the Hartree shift in more detail. In the future, it will be 
interesting to consider approximative methods that combine the two, the BCS theory 
and the Hartree shift. 

Solving the k frequency integral in Eq. ^ (by closing the integration path in 
a semi-circle around the upper half of the complex plane) yields contributions from 
poles in the Green's function G but also from possible poles in T that are located in 
the upper half of the complex plane. The self-energy can thus be separated as 

E(P) = SHartrco + Spairs, (4) 

where the Hartree part Enartree of the self-energy corresponds to the part given by 
the pole in the Green's function and the pair part £ pa irs is given by the poles in the 
T-matrix. This is the definition of the Hartree self-energy used in this work: it is the 
part of the self-energy arising from the pole in the Green's function in Eq. In 
addition, by the Hartree shift I will refer to the real part of the Hartree self-energy 
Re Snartree- Mathematically speaking this separation to Hartree and pair self-energies 
is very easy to do as the contributions from the poles add up exactly in such a way 
when doing contour integration. Physically speaking the poles in the Green's function 
give the occupation numbers and, eventually, the atom density. The poles of the T- 
matrix correspond to many-body bound states or molecular states if the two-body 
interaction potential supports such states. The simplest theory for describing these 
many-body bound states would be the BCS theory, in which case the Cooper pair 
formation is signaled by the appearance of a pole in the many-body T-matrix at zero 
center-of-mass momentum and energy. However, the pole of the Green's function 
is neglected in the BCS self-energy for strong interactions, as the dominant term is 
considered to be the pole in the T-matrix describing the Cooper pair formation. The 
above definition for the Hartree self-energy is very similar to the self-energy used in 
Ref. [22], and notice also that it yields Eq. ^ in the weakly interacting limit when 
the T-matrix is replaced by the constant Antra/ni. 

3. Many-body scattering T-matrix at unitarity 

I will start by reviewing the derivation of the many-body scattering T-matrix. While 
the final result is already well known and widely used, the actual derivation is not 
often shown, with the calculation often starting from the two-body scattering T-matrix 
instead of the actual atom-atom interaction potential. However, going through the 
actual derivation may help to understand the physics of the Hartree self-energy and 
where the various features in the many-body scattering T-matrix arise from. The 
detailed calculation can be found in |Appcndix A[ but here in the main text I will do 
the key approximations and show only the main results. 
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3.1. Inelastic scattering amplitude 

The first key result from the scattering calculation in jAppcndix A| is the elastic 
scattering amplitude for a short range potential 

/(s,k) = /(fc)« — ioxk = s (5) 
1 + ika 

which describes the scattering of two particles with the incoming and outgoing 
(relative) momenta k and s, respectively. For short range potentials, the scattering 
length a is enough to parametrize the two-body scatterings in vacuum. As described 
in |Appcndix A| I will assume that the inelastic scattering amplitude has the same 
form, i.e. /(s, k) = f{k) for all s,k. While such assumption is not true in general, it 



is consistent with the generalized optical theorem for inelastic scatterings in Eq. (A.9) 
at unitarity a — > oo. 

At low energies, k sa 0, the scattering amplitude is constant —a, but the 
momentum dependence and the associated decay of the real part of the scattering 
amplitude, 

implies the well known fact that the interaction strength decreases when the relative 
momentum of the two particles increases. For interaction potentials with a finite range, 
the decay is even faster. This means that a very fast atom does not feel interaction 
effects and all kinds of pair formation and Hartree energy shifts are suppressed. The 
idea of a constant Hartree shift thus holds only for a limited range of momenta, with 
the characteristic scale given by the inverse of the scattering length I /a. Obviously, 
in strongly interacting gases the condition ka <C 1 is easily breached. 

The momentum dependence of the Hartree energy shift has been observed in 
Bragg spectra of strongly interacting Bose-Einstein condensates [5^1 E01 EI] - However, 
compared to strongly interacting fermionic gases, the associated momenta k and 
scattering lengths a in that experiment were still rather meager. In this work I consider 
only the two limits: the weakly interacting limit /(s,k) ~ —a and the unitarity 
/(s,k) = i/k, but the theory is quite easily generalized to intermediate interactions 
as well. 

3.2. Two-body scattering T-matrix 

The scattering amplitude describes stationary scattering states but in order to 



calculate the self-energy the two-body scattering T-matrix in Eq. (A. 33 1 is needed 



4ttH 2 

r (s,k;F)= /(s,k) (7) 

m 

4Trh 2 \ 2 f dq ( 1 



V m ) J {211)* \2p - 2e p - 2e q 2e q ~ 2e k - ir) J 

The physical interpretation of the two-body scattering T-matrix is very simple. 
The scattering amplitude /(s,k) obtained from the standard scattering calculation 
describes stationary solutions of the scattering problem, and thus the energies of 
the incoming and outgoing particles are equal (there is only one 'energy' for each 
eigenstate of the Hamiltonian) . This implies that the energy is not a free parameter 
in the scattering calculation. The two-body scattering T-matrix To now generalizes 
this scattering amplitude into the case in which the energy is a free parameter. This is 
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important when considering for example scattering processes where the particles can 
exchange energy with the environment (say, a probing radio-frequency or laser field) 
but also when the energy is not a well defined concept such as in the case of virtual 
intermediate states. 

In the weakly interacting limit the two-body scattering T-matrix simplifies into 

r (s,k;P)= - 4 ^+0(a 2 ). (8) 
m 

At unitarity it becomes 

4irh 2 i (4irh 2 \ 2 f dq 1 / 1 1 \ 

r (s,k;P) - ~ — — j J —- y——— + ———J . (9 ) 

The first important point to notice is that the two-body scattering T-matrix does not 
depend on the momentum s at all, i.e. the approximation done for the scattering 
amplitude, /(s, k) ss f(k), carries on to the two-body scattering T-matrix. 

The integrand is spherically symmetric and the integral is easily evaluated 

r (s,k;P) = 




where in the last line I have done the analytic continuation po — e + in with an 
infinitesimally small imaginary part r\ — > 0. The two-body scattering T-matrix is thus 
either real or pure imaginary, depending on the sign of e p — e. Notice also that the 
two-body scattering T-matrix does not depend on either momenta s or k but only on 
the center-of-mass momentum p and energy pq. However, one should notice that this 
property is somewhat illusionary: for an elastic scattering the two-body scattering 
T-matrix is identical to the scattering amplitude f(k), which depends on the relative 
momentum k instead of the center of mass momentum p. Indeed, for an elastic 
scattering e — e p — and the two-body scattering T-matrix regains the form used for 
the scattering amplitude. 

3.3. Many-body scattering T-matrix 

The many-body scattering T-matrix can be expressed using the two-body scattering 



T-matrix as in Eq. ( |A.39| ) 

r(s,k;P) = r (s,k;P)+ j ^r (s,q;P)x(q;P)r(q,k;P). (11) 

Using the knowledge that the two-body scattering T-matrix depends only on the 
center-of-mass four-momentum P allows writing this as 

r(p) = ^5 , (12) 

i-r (P)/^x(q;P)' 



Combining this with the result obtained for Tq in Eq. ( 10 1 yields 



nP) = -m (13) 
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where 



where V signifies the principal value integration. This is exactly the same as the 



g(P) = i t^G(P + Q)G(P -Q)+ — ^ — , (14) 

(27r) 4 J (27r) J 2p - 2e q - 2e p 



standard T-matrix at unitarity a — > oo jy] , and as pointed out at the end of Appendix 
[K\ the many-body scattering T-matrix yields the standard result whenever the two- 
body scattering T-matrix does not explicitly depend on the relative momenta s,k. 

In order to proceed further, I will make an ansatz for the Green's function G(P). 
The choice of the ansatz is probably the most important approximation in this work 
and hence it must be chosen with care. I assume that the Green's functions have the 
form 

G(P)= . ni£p \, (] + l ~ n{e iy M 
Po - e p - h \P) Po-fp- h(p) 

where n(e p ) is the occupation probability for a single-particle state with energy e p 
(or momentum p). The first term, the advanced Green's function, describes a hole 
excitation with energy (dispersion) e p + h*(p) and the second term, the retarded 
Green's function, describes a particle excitation. The complex conjugation of h(p) in 
the hole branch guarantees that the poles corresponding to holes remain in the upper 
half of the complex plane and the particle poles in the lower half (the imaginary part 
of h(p) must thus be negative). The key approximations in this ansatz are that the 
self-energies h(p) for the advanced and retarded parts are the same and that they do 
not depend on the frequency pq. These approximations exclude a lot of physics from 
entering the formalism, perhaps the most important piece missing being the formation 
of an excitation gap. Indeed, the BCS Green's function would require a self-energy of 
the form S(P) = — Go(— P)A 2 , where A is the BCS excitation gap and Gq(P) is the 
non-interacting Green's function, which clearly depends on the frequency po. However, 
the approximations fit with the picture of a momentum-dependent Hartree shift, with 
each Green's function pole being simply shifted by the amount Reh(p). By allowing 
complex values of h(p) the formalism also incorporates lifetime effects, broadening 



the spectral function. The ansatz in Eq. (15) is thus sufficient for describing Hartree- 



like effects. Taking into account also the pairing effects would generalize the present 
theory in the direction of the pair-fluctuation theory in Ref. [15] . However, here I will 
concentrate on the Hartree effect and ignore the pairing contribution. 

The ideal zero-temperature Green's function can be obtained from the ansatz by 
choosing n(x) = 1 — 9(x) and non-interacting finite temperature Green's function by 
using the Fermi-Dirac distribution. Below I will also consider occupation numbers of 
the form given by the Bogoliubov coefficients 

n(e P ) = vl = if 1 - ^ ) (16) 



in order to have an ansatz that satisfies the Bethe-Peierls condition for the two-body 
physics jTUJ [33 - However, despite the similarity with the BCS Green's function, this 
does not result in the emergence of an excitation gap. Instead, the parameter A is 
simply a measure of the broadening of the Fermi surface. In order to emphasize this 
difference with the true BCS Green's function, but also to show the purpose of the 



ansatz, I will below refer to the ansatz in Eq. (16) by the 'BPC ansatz'. The BPC 
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Figure 1. The real and imaginary parts of T(s, d) as a function of the center-of- 
mass momentum s and the relative momentum d calculated for a zero-temperature 
Fermi sphere. 



ansatz will be discussed further in Section [6j As will be seen later, the Hartree energy 
shifts are not very sensitive to the exact form of the occupation number ansatz. 

The regularized pair susceptibility can now be evaluated using contour integrals 
yielding 



9( p) = / 



n(-e+)n(-e_) n(e+)n(e_) 



(17) 

. 2p -£+-£_ 2pg - e* + - e*_ 2p Q - e + - e_ J ' 

where e± — e p ± q and e± = e± + /i p ± q . This, together with Eq. (13) giving 
the connection to the many-body scattering T-matrix, is the key result needed for 
calculating the Hartree self-energy. The appearance of the dressed Green's function 
G{P) in the definition of the many-body scattering T-matrix implies an iterative 
solution. Indeed, the present theory is fully self-consistent in the sense that only 
dressed Green's functions are used both in the T-matrix and in the self-energy. 
Fortunately, due to the simplicity of the Green's function ansatz and because of a 
rapid convergence of the iteration, the resulting numerical task is very easy to solve. 
As will be described below, I will also consider an alternative scheme in which the 
imaginary part of the self-energy is neglected. 

3.4- Numerical solution of the many-body T-matrix 

In order to obtain a better physical picture of how the many-body scattering T-matrix 



looks like, I will for the moment neglect the Hartree self-energies h{p) in Eq. (17). As 



will be seen below, for the determination of the Hartree self-energy, it is sufficient to 
consider elastic scatterings for which P = ( , £p ^ £fc ) . The many-body T-matrix 
becomes now a function of only two parameters, the magnitudes of the center-of-mass 
momentum s and the relative momentum d. 
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Figure 2. The real and imaginary parts of F(s,d) as a function of the center- 
of-mass momentum s and the relative momentum d calculated for a finite 
temperature T = 0.5 Tp Fermi sphere. 




Figure 3. The real and imaginary parts of d) as a function of the center- 
of-mass momentum s and relative momentum d calculated for the BPC ansatz 



occupation probabilities n(x) = J ( 1 . x B == ) with fi = 0.59 Ep and 

\ vC^ - m) 2 +A 2 

A = 0.69 E F . 
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Figure 4. The real part of F(k z ,k r ) as a function of the atom momentum 
k = (k z ,k r ) for an atom with momenta p = (top-left), p = 0.50 kp (top- 
right), p = 1.0 kp (bottom-left), and p = 2.0 kp (bottom-right). Shown is the 
absolute value of the real part, calculated assuming the zero-temperature Fermi 
distribution for the occupation numbers. 



Fig. [T] shows the real and imaginary parts of the many-body T-matrix T(s,d). 
The most striking feature is the divergence of the real part for small relative momenta 
d at sufficiently high center of mass momenta s > kp. This is caused by the 1/d 
divergence of the two-body scattering T-matrix Tq but one should remember that Tq 
was a purely imaginary quantity for elastic scatterings, whereas also the real part of T 
diverges. The imaginary part of T reproduces this behavior too, but one should also 
notice the suppression of the divergence for low center-of-mass momenta s. There is 
also a peculiar feature appearing around d = 0.849 kp, s = 0.546 kp. The physics of 
this feature will become clear shortly. 

Fig. [2] shows the many-body T-matrix for a high temperature T = 0.5 Tp Fermi 
sphere, showing the broadening of the features when the Fermi surface becomes less 
well defined. In additio n, F ig. [3] shows the many-body T-matrix calculated assuming 

the BPC ansatz in Eq. (16 1 for the occupation numbers n(x) — k ( 1 



l 

with fi, = 0.59 Ep and A = 0.69 Ep- Because of the similarity between the results 
obtained using this ansatz and the high temperature ansatz I will below concentrate 
only on two cases: the zero temperature (with a well defined Fermi surface) and the 
BPC ansatz. 
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Figure 5. The real part of F(k z ,k r ) as a function of the atom momentum 
k = (k z ,k r ) for an atom with momenta p = (top-left), p = 0.50 /cf (top-right), 
p = 1.0 fcp (bottom- left), and p = 2.0 fep (bottom-right). Shown is the absolute 
value of the real part, calculated assuming the BPC ansatz for the occupation 
numbers. 



Despite the simplicity of the T(s, d) plots, it is still rather hard to see from these 
plots which interaction channels are actually dominating. Fig. [4] shows the many-body 
scattering T-matrix as a function of actual atom momenta k for a test atom with a 
momentum p, shown for different values of p. The strong peak observed for an atom 
with high momentum p is easily understood as arising from the 1/d divergence in the 
high center-of-mass momentum scatterings. In contrast, the scattering channels for a 
low momentum atom (p = Okp) show a very distinct ring with a radius of roughly 
k = y/2kp. This ring corresponds to the small feature seen at the low center-of- 
mass momentum in Fig. [T] The radius of this ring, and the corresponding remnants 
of the ring for intermediate momenta (p = 0.5 kp), suggest that the atoms would 
interact most strongly with such atoms that allow on-shell scatterings onto the Fermi 
surface. Indeed, two atoms with momenta p = and k = y/2kp can scatter right on 
to the Fermi surface with magnitudes of momenta q — fcp. Likewise, an atom with 
momentum p — 0.50 fcp has a strong interaction peak with an atom with momentum 
k m (0.805, 1.05) kp. The total kinetic energy of these two atoms equals 2 Ep, allowing 
them to scatter on to the Fermi surface. However, interference effects from different 
scattering channels play an important role resulting in the appearance of clear maxima 
and minima on and inside the ring. One should notice, however, that the many-body 
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scattering T-matrix simply tells how strongly the atoms with momenta p,k would 
interact regardless of whether the two atoms exist or not. Indeed, for an atom inside 
the Fermi sphere, the strongest interaction channel is always with a momentum state 
outside the Fermi sphere. 

Broadening the Fermi surface by either a finite temperature or by the BPC ansatz 
makes these features smoother as seen in Fig. [5] For low momenta p, a clear ring is still 
seen, with the radius k ss 0.77 fcp corresponding to the kinetic energy h 2 k 2 /2m « 2/x. 
The imaginary parts of the many-body scattering T-matrix look qualitatively similar 
to the real parts and thus they are not plotted here. 



4. Hartree self-energy 

The main quantity of interest in this work is the Hartree self-energy. As discussed 
in Section [3] this requires a self-consistent iterative solution and hence the inclusion 



of Hartree self-energies in the many-body scattering T-matrix in Eq. (17 1. The self- 
energy was defined in Section [2] as 



E(P) = -t J -^G{K)T{K,P-P,K) 



dK 

G(k,fc )r(d,d;5), (18) 

(2tt) 4 

where S = (s, p "+ k ° ^ ; d = ~(p — k), and s = ~(p + k). Doing the k frequency 

integral yields contributions from the pole of the Green's function G but also from 
possible poles in the upper half of the complex plane in T. The first part was defined 
as the Hartree self-energy. Since we are interested in the energy of a particle with 
momentum p, we will evaluate the self-energy at frequency po = e p . Notice that this 
results in an implicit equation for the Hartree self-energy, as the point where the value 
of the self-energy is to be evaluated depends on the value of the self-energy through 
Sp 6p hp. 

We can now write the Hartree self-energy as 

„ (e|l )r(d,d;B,— 



Uari,'ee(P,Ep) = / n(e k )T (d, d; &, — ) (19) 



Notice that only the particle branch of the Green's function G(K) is relevant as the 
pole corresponding to the hole branch is on the wrong side of the complex plane. 
The Hartree self-energy is a complex quantity with the real part describing the actual 
Hartree energy shift and the inverse of the imaginary part yielding the lifetime of the 
corresponding excitation (hole or particle excitation). However, as will be seen below, 
omitting the pairing part of the self-energy causes problems with the imaginary part 
of the total self-energy namely that the imaginary part of SHartrco is always negative. 
This means that the poles corresponding to both particle and hole excitations are 
shifted upwards in the complex plane in violation of the Lehmann representation 
which states that the imaginary part of the self-energy must change sign at the Fermi 
momentum [33] . 

This problem however does not manifest in the present model because of the choice 



of the Green's function ansatz in Eq. ( 15 ) where the imaginary part of the self-energy 
in the hole and particle branches have opposite signs by design. The correct behavior 
of the poles is thus forced, and since the hole and particle branches also use the same 
self-energy the problem does not appear. Being aware of these problems with the 
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imaginary parts of the self-energy, I will also consider an alternative path by relaxing 
the requirement of self-consistency and omit the imaginary parts of the self-energy. 
Such approach certainly underestimates the effect of lifetime broadening (by neglecting 
it altogether) while the fully self-consistent approach most likely overestimates it since 
the imaginary part of the self-energy should vanish at the Fermi surface [33] but it 
does not in the present model. We will use the self-consistent approach in most of the 
results below because the non-zero imaginary part of the self-energy broadens spectral 
peaks making numerical integrations easier. Both approaches give qualitatively similar 
results and they can probably be considered as the upper and lower bounds for the 
Hartree energy shifts and the excitation lifetimes. 



4-1- Weakly interacting limit 



The problems with the imaginary part of the self-energy can be easily observed by 
first considering the weakly interacting (and low momentum) limit. This calculation 
is effectively a review of the Galitskii's calculation for the lifetime and energy shift of 
quasiparticle excitations. 

For weakly interacting atoms the scattering amplitude can be approximated as 

/(s, k) =w —a + ika 2 , 

yielding the two-body scattering T-matrix 



r (s,k;P) 



-ka 2 



Anh 2 a ( Airtfa 
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m \ m 



/Airh 2 a 
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dq 



dq 
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1 



1 



2p - 2e p - 2e q 2e„ - 2e k 



(2tt) 3 \2 Po ~2e p -2e q 2e q - 2e k 



where V signifies the principal value integration and terms of order 0(a 3 ) have been 
dropped out. Notice that this does not depend on s but it does depend on k. 

Assuming ideal zero-temperature Green's functions G(K) = Go(K) for the 
intermediate scattering states, the many-body scattering T-matrix becomes (here too 
omitting terms 0(a 3 )) 

dq 



r(s,k;P) = r (s,k;P 



(2tt) 3 
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(20) 



(21) 



n(-e p+q )n(-e p _ q ) n(e p+q )n(e p _ q ) N 



2p - 2e p - 2e q 



and the Hartree self-energy is now 



-"Hartree 



dk 
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n(e fc )r d,d;s 
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dq 
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'fl (|s + q|-l)fl(|s-q|-l) _ fl(l-| a + q|)0(l-|s-q|) V_ 

d 2 — q 2 + irf d 2 — q 2 — irf d 2 — q 2 t 

where r( — mr//h 2 and the relation n(e/.) = 0(Ep — e/.) = 9(1 — k/k-p) has been used in 
order to allow a direct comparison with the results in Ref. |33j . The dominant (first) 
term is the standard weakly interacting limit Hartree shift and the rest of the terms 
are corrections to this weak interaction result. However, a careful comparison with 
Ref. [33 shows that the result is missing terms proportional to 8(k — 1) and it has 
an extra term proportional to 8(1 — k). The reason for the discrepancy is the neglect 
of the pole in the T-matrix when doing the k$ integration. However, as discussed in 
Section [2] such effects do not enter the Hartree part of the self-energy but instead 
would go in to the pair part S pa i rs . Notice, in particular, that the imaginary part of 



the Hartree self-energy in Eq. ( 23 1 does not change sign at the Fermi surface and hence 



the imaginary part of the Hartree self-energy for hole and particle excitations would 



have the same sign. The Green's function ansatz in Eq. ( 15 ) fixes this problem by 
demanding that the self-energy for the hole branch is obtained by taking the complex 
conjugate of the particle branch self-energy. 



4.1.1. Unitarity The Hartree self-energy at unitarity can be calculated from Eq. ( 19 1 

dk 



J Hartrcc 



(P) 



(2tt) 3 
dk 

W) 



n(e k )T [ d,d;s, 
3 n(e k )T 



(23) 



It is instructive to consider first the case in which the many-body corrections are 
ignored. This can be achieved by approximating r ~ Tq and e p ~ e p , yielding 



Hartree 




Hp+ k ) 



(24) 

The lesson here is that the Hartree self-energy does indeed involve the relative 
momentum p — k instead of the center-of-mass momentum. This is a result of the 
way the center of mass quantities are formed, namely that the energy difference 
€p ~^ £k — e s = e^. The same energy difference enters also the many-body T-matrix 
and thus also it will depend on the momentum difference d. Notice that the Hartree 
self-energy obtained using the two-body scattering T-matrix Tq is purely imaginary 
and thus the atoms would not experience any energy shifts but excitations would only 
have a finite lifetime. The next step is to evaluate the Hartree self-energy by including 
the many-body corrections. 



4-1.2. Numerical solution of Hartree self-energy at unitarity Fig.[6]shows the Hartree 
self-energy Euartree(-P) as a function of the atom momentum p (and for energy po = e p ) 
for both the BPC ansatz occupation distribution and the zero-temperature Fermi 
distribution. The Hartree shift is strongly momentum dependent, although for low 
momenta p < fep the dependence is not very strong. However, beyond the Fermi 
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Figure 6. The real and imaginary parts of the Hartree self-energy £Hartree(p> £p) 
as a function of the atom momentum p for the ideal Fermi sphere and the BPC 
ansatz. 

momentum, the Hartree shift exhibits a 1/p decay. It is also worth noticing that 
the energy shift is higher at low momenta for the BPC ansatz than for the well- 
formed Fermi sphere. This is because of the appearance of holes in the Fermi sphere, 
relaxing the effect of the Pauli blocking. This effect can be seen also by increasing the 
temperature in the Fermi distribution. 

Shown is also the imaginary part of the Hartree self-energy with a pronounced 
peak slightly above the Fermi momentum. As discussed in Section[4j this behavior is at 
odds with the Galitskii's result for a weakly interacting gas, where the imaginary part 
vanishes at the Fermi momentum. On the other hand, the pair part of the self-energy 
E may also provide an imaginary part that can cancel at least partly the imaginary 
part in the Hartree self-energy Snartree- 

Fig. [7] shows the real and imaginary parts of the Hartree self-energy calculated 
using the alternative scheme in which the imaginary part of the self-energy is omitted 
altogether. Notice that even though the imaginary part of the self-energy is not 
included in the iteration, its value can still be calculated and hence shown in the figure. 
While the numerical values for the two schemes do differ somewhat, the qualitative 
behavior is similar for the two approaches. As argued in Section [4] the 'true' behavior 
is expected to be somewhere between these two limits. 

5. Hartree self-energy in a spin-imbalanced gas 

The present method can be easily generalized for studying polarized, or spin- 
imbalanced, gases. This can be accomplished by having spin-dependent distribution 




Figure 7. The real and imaginary parts of the Hartree self-energy £Hartree(p, £p) 
as a function of the atom momentum p for the ideal Fermi sphere and the BPC 
ansatz, calculated using the semi-self-consistent scheme in which the imaginary 
part of the self-energy is neglected. 



functions n a {x) and Hartree shifts h a . p , with a G {t?i}- The only equation that needs 
to be modified is the g(P) function in Eq. (171, which becomes 



g(P) = 



dq 



«t(- 



tp+q 



)%(- £ p-q) «t( e P+q)«i( e P r q) 



2po 



where e a .± — e± + h a _±. Notice that g{P) is symmetric with respect to the change 
'fi — >l- and the functions g and T are independent of the spin. 



2p5 



1 

2e - 2e„ 



2e 



,(25) 



The Hartree self-energy Eq. ( 19 ) for atoms with the spin f becomes 



n T (e fc )r d, d;s 



(26) 



(27) 



and for spin J, atoms it becomes 

Hartree ~ j 

Fig. [8] shows the Hartree self-energy for a spin-imbalanced gas with the 
polarization P w 0.5 at the zero temperature. What is striking is the increase of 
the Hartree energy shift experienced by the minority component as compared to the 
balanced gas (the number of majority component atoms is the same for balanced and 
spin-imbalanced cases). This is the result of the relaxation of the Pauli blocking, 
as the low momentum minority component atoms have access also to low momentum 
scatterings with majority component atoms close to the Fermi surface. In contrast, in a 
balanced gas such low momentum scatterings are blocked and the 'lowest momentum' 
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Figure 8. The Hartree self-energy of a polarized (P = 0.5) zero-temperature 
Fermi gas. The real part is shown in the left and the imaginary part on the right. 
The Hartree self-energy for the balanced gas is shown for reference. 



scattering would amount to lifting the atom onto the surface of the Fermi sphere. 
Obviously the corresponding energy shift is decreased for the majority component as 
there are fewer minority component atoms to interact with. 



5.1. Polaron limit 

An interesting special case is the extreme polarization limit P — > 1 where one has 
only a single minority component atom as an impurity interacting with an ideal Fermi 
sphere of majority component atoms. Such a setting is predicted to create polarons 
and it has been already widely studied both theoretically [351 [Ml [37J [Ml ESI ED] and 
experimentally [HI H3J H2] . In the polaron limit one would expect that the pair 
part of the self-energy has the least effect and the prediction from the Hartree part 
would be thus of particular interest. 

However, first I would like to make a connection to the polaron theory in Rcf. [37 , 
which also led to the same analytical formula as the variational Chevy ansatz in |36) . 
In Ref. [37] the self-consistent iteration was cut after the first iteration step, resulting 
in a theory where the self-energy of a particle with momentum p was evaluated at 
energy e p + h^ p , but the Hartree shift h^, p was not fed back into the many-body 
scattering T-matrix, yielding as the polaron energy at unitarity a — > oo 

Notice the implicit form of the solution, which is the result of the important assumption 



Hartree shift in unitary Fermi gases 18 




'■ 9 0.5 1 1.5 2 2.5 3 "°- 7 0.5 1 1.5 2 2.5 3 

Momentum p/k F Momentum p/k F 

Figure 9. The real and imaginary part of the Hartree self-energy of a single 
minority component atom immersed in a zero-temperature Fermi sea of majority 
component atoms. Shown are results from both the fully self-consistent scheme 
and from the semi-self-consistent scheme in which the imaginary parts of the 
self-energy are neglected in the iteration. 



that although the Hartree shift is not fed into the many-body scattering T-matrix, it 
does enter in the energy E at which the energy shift is evaluated at. In the language of 
this work, it corresponds to using bare Green's functions Go in Eq. ( 14 ) but evaluating 
the Hartree energy shift in Eq. ( 19 1 at a shifted energy e p = e p + h p . While this appears 
like a strong assumption, the results (with the zero-momentum polaron energy being 
roughly — 0.606 Ep) were in very good agreement with Monte Carlo results [35]. The 
Hartree self-energy described in this work simplifies into Eq. ( 28 1 when doing the same 
approximations but it also allows a more rigorous analysis of the justification and the 
meaning of these approximations. 

The simplest generalization from Eq. (28) would be to include the Hartree shift 
in the many-body T-matrix but assume that it is constant, i.e. it does not depend 
on the momentum h±p = h±. This corresponds to adding the Hartree shift to the 
single-particle dispersions of the spin-| atoms in the second term in Eq. (28 1, yielding 



E 




-p+q— k 



h i 



E 



(29) 



Notice that only one of the energies in the second term describe a spin-^ atom and 
hence the shift h± appears in the denominator only once (the energy shift of majority 
component atoms is assumed to be zero). However, the second term in the right is 
the Hartree energy shift hi , and one obtains E = e.„ + h^ 1 which upon substitution to 
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Eq. p9|) yields 

„ / „ i i \ 

(30) 




This corresponds to the approximation in which the Hartree self-energy in Eq. ( 19 ) is 



evaluated at energy e p . Solving this equation is actually easier than Eq. (281, yielding 
as the zero-momentum polaron energy E « — 0.917 Ep. The excellent agreement with 
the Monte Carlo result is thus lost. 

Considering these issues from the point of view of the scattering theory presented 
in this work, the reason why the Chevy ansatz works so well is that the self- 
energy is strongly momentum dependent, so that most of the intermediate scattering 
states have vanishing self-energy. Indeed, due to the short range of the interactions 
and the corresponding minor importance of low energy states in the scatterings, most 
of the intermediate scattering states are only weakly interacting with the majority 
component Fermi sphere. This in turn limits the level of correlations needed for an 
efficient description of the system and the simple model provided by the Chevy ansatz 
works well. Contrast this with a system in which the intermediate scattering states 
are limited by a low momentum cutoff (for example due to a long effective range of the 
interaction potential or due to a band gap in the single-particle excitation spectrum 
in a deep optical lattice). In such case approximating these intermediate states with 
a non-interacting Green's function would not work as well. 

Solving the self-energy in a self-consistent manner by doing the required Hartree 
shift iteration yields the Hartree shift shown in Fig. [9j It shows the polaron self-energy 
at unitarity as a function of the polaron momentum p and, as expected, it exhibits a 
strong momentum dependence with a 1/p decay at high momenta. At low momenta 
p < kp, the energy shift is roughly constant with magnitude approximately —0.61 Ep, 
which is in good agreement with predictions from other theories [35l [36l [37j [3Sl S3 [39] 
and experiments [?! ] l42 j l43l 144] . 

From the momentum dependence of the self-energies one can also obtain the 
effective mass to* for the polaron [37] 

as 



m*=m - j - . (31) 

1 + <9(p 2 /2m) 

Ignoring the possible frequency dependence of the self-energy, I obtain the effective 
mass to* rj 1.04 m at unitarity, which agrees well with Monte Carlo results [35l [45] 
and the analysis of an experiment in Ref. [41] . although the value differs somewhat 
from the predictions in Refs. [37J [35] and observations in new experiments [121 WQ . 
Neglecting the imaginary parts of the self-energy gives the zero-momentum polaron 
energy — 0.70 and the effective mass to* 1.47 to. 

6. Contact 

The self-energy yields also the momentum distribution through the relation [33] 

n k = J^Gfako), (32) 

but for this one would need to determine the Green's function and hence also the 
self-energy for all frequencies ko- However, I propose that for a dilute gas the 
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high momentum asymptote of the distribution can be calculated using a simple 
perturbation theory. 

Consider an adiabatic switching on of the interactions, H(t) = Hq + Ve~ vt for 
t < and 77 — > 0, starting from the ground state of Hq. According to Gell-Mann and 
Low this yields the ground state of H as 

|V) =S(0, -00) I Vo>, (33) 

where \4>q) is the ground state of Hq and the S-matrix 

S(t,t')=Te- i ti dtlV ( tl \ (34) 

where V(ti) is the interaction part of the total Hamiltonian expressed in the interaction 
picture. 

The expectation value of the occupation number for the spin state f i n some 
momentum state n q/ j- = (fl , \c cl *Cq ! f\ip) in the interacting ground state is now 

«q,t = |Cq )t S'(0, — OO)!^)! • (35) 

Motivated by the arguments of Shina Tan [47l 08], I assume that only pair- wise 
correlations are relevant in the high momentum limit q 3> kp , and thus the observation 
of a particle at the momentum q implies that some pair of atoms (with momenta p and 
k) in the initial state |^>o) have scattered to the states q and p + k q. The probability 
amplitude for such scattering is provided by the many-body scattering T-matrix T(P) 
(which, according to Section [3] depends only on the center-of-mass four-momentum 
P) divided by the energy change due to the scattering, i.e. the probability is 



P t (k,p;q) = 



2 - 2 



n k n v 



r(p) 

2e Q 



n k n p , (36) 



£ p+k-q + £ q ~ £ P ~ £ k 

where the occupation numbers n k and n p give the probability that the momentum 
states k and p were initially occupied and I have assumed that q 3> p,k and that 
e q k, e q (based on the observation that the Hartree energy shift vanishes for large 
momenta q). This is similar to the first-order perturbation theory result for the 
transition probability assuming the coupling strength T(P) and the energy difference 
2e q . However, notice that since the coupling strength is in reality a T-matrix, it 
includes scatterings to arbitrary order but the correlations are limited to the two- 
particle level. 

However, the observation of a particle at the momentum q tells only that some 
atoms have scattered but does not specify the original momenta of the atoms p, k. 
On the other hand, since higher order correlations are ignored, there are holes left at 
the original momenta p, k. Another trajectory leading into observation of a particle 
at momentum q would be one in which two atoms from some other momentum 
states p',k' scatter to states q, p'+k' — q. However, the two trajectories yield a 
zero overlap when calculating the expectation value since the final states are mutually 
orthogonal due to the different configuration of holes. Hence the two trajectories do 
not interfere coherently and one expects that the different trajectories need to be 
summed incoherently when calculating the expectation value of the observable c q c q . 
Therefore I assume that the occupation number at some high momentum state q 3> kp 
is given by 

n i.t ra 772 ZJ^'P^I ~ ^^72 > ( 37 ) 

k,p 1 
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The contact [47j [48], defined as the limit C = lim^oo q n, 

C =(S) 2 ^^ |r(P)|2 ^' fe ^- (38) 
k, P 

The value of the contact depends strongly on the low momentum occupation numbers 
rik both directly but also indirectly through the many-body scattering T-matrix F. 
Notice that this equation differs from the one used for example in Refs. [321 US] where 
the occupation numbers were defined through self-energies. It is instructive to consider 
the weakly interacting limit first. Approximating V pa 4 yields for the contact 



/m\ 2 1 / 4irh 2 a x 



Cw liJ V^(^r) ^,kH, P = ^(nar, (39) 

^ ' k,p 

where n is the total atom density of a balanced gas. This is the exact result for a 
weakly interacting gas [3H1 EE] • 

Another simple system for which one can easily determine the contact is the 
polaron system where the majority component atoms form an ideal Fermi sphere 
n^,k = 0(kp — fc) and the minority component atoms consist only of a single impurity 
at the zero momentum n^ p — S p q. The contact becomes now 

ir(jj,s*±a*)i%,. (40) 

k 

Notice that this quantity depends on the Hartree shifts at all momenta because these 
affect the intermediate states in the many-body scattering T-matrix. At unitarity I 
obtain C — 0.10 fc-^pfcj F (or C — 0.13 fc^pArj' F when neglecting the imaginary part 
in the self-energy). This is slightly higher than the value C pa 0.08 fct,pA;| F obtained 
in [ID]. I am not aware of experimental results for the polaron contact. 

The contact for a balanced system has been determined experimentally using 
various methods but the problem with the present theory is that it does not give 
access to occupation numbers for low momenta k ~ fcp. The simplest assumption 
is that both atomic components form an ideal Fermi sphere with occupation numbers 
rife = 9(kp — k). Such assumption gives for the value of the contact C pa 0.064 /cp. 
Omitting the imaginary parts in the self-energy yields C ~ 0.14 fcp. Obviously, any 
finite value for the contact is incompatible with the assumption of a Fermi-Dirac-likc 
distribution for the occupation numbers as it has an exponentially decaying tail in 
violation with the Bethe-Peierls condition. The assumption can be easily improved 
by making an ansatz for the occupation numbers using a functional form that has the 
correct asymptote. For this end I use the BPC ansatz 

where fi and A are considered fitting parameters with fi chosen so that the atom 
numbers are fixed and A is obtained from the value of the contact C = limfc^oo q 4 n q = 
(|rj) 2 A 2 . This requires an iterative solution so that the contact obtained from 
Eq. ( [38] ) matches the contact from Eq. (41). Using this ansatz I obtain A sw 0.49 Ep, 



/i w 0.79 E-p and C ps 0.060 fcp. In addition, the Hartree energy shift at the Fermi 
momentum equals — 0.34£p. Neglecting the imaginary parts of the self-energy yields 
A ss 0.64 Ep, fi « 0.66 Ep, C pa 0.103 fcp, and the Hartree shift at the Fermi 
momentum —0.46 -Bp. 
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The difference in the vaiues from the two approaches for the self-consistency (as 
discussed in Section [4]) can be understood in terms of the hole lifetime. Including the 
imaginary parts of the self-energy in the Green's functions of the scattering particles p 
and k will decrease the expectation value of the scattering probability since the system 
cannot return to the initial state if the holes have decayed. In the extreme limit of a 
very short lifetime, or very large imaginary parts of the self-energy, all scatterings will 
be blocked analogously to the quantum Zeno effect. 

Considering the simplicity of the present method, these figures are in surprisingly 
good agreement with Monte Carlo results [50] (yielding C — 0.1147(3)fcp), self- 
consistent Green's function method [T7] (yielding the pairing gap A = 0.46 Ep, the 
Hartree shift U = —0.50 E-p, and the contact C — 0.102 k F ), and experiments (Ref. [51] 
yielding C 0.12 k F and Ref. [52] yielding C « 0.096 k F ). 



7. Discussion on the model 



The key approximations in this work are the Green's function ansatz in Eq. ( |15| 
and the neglect of the pairing part of the self-energy as described in Section [2] The 
justification of both of these approximations is based on the argument that they are 
sufficient for describing the Hartree energy shift but one should not make too far 
reaching conclusions based on this model. However, as seen by comparing the results 
obtained for the zero-temperature Fermi sphere ansatz and for the BPC ansatz in 
Figs. [6] and [7J and also the values obtained for the contact in Section [6j the values 
of the Hartree self-energy do not have a very dramatic dependence on the ansatz. A 
distinct Fermi surface does make features sharper, see Fig. [7J but qualitatively the 
behavior is similar and the actual numerical values close to the Fermi momentum 
are in the same range. Thus one would expect that changing the actual form of the 
BPC ansatz would not change the values much. One interesting possibility would 
be an ansatz that combines a step-like behavior at the Fermi surface (albeit with 
a step smaller than unity) and a proper 1/k high momentum asymptote or, even 
better, describing the low momentum regime with the finite temperature Fermi-Dirac 
distribution and continuously merging this with a 1/fc 4 tail. The result of such an 
ansatz would be a Hartree shift lying somewhere between the two ansatzes used in 
this work. 



The BPC ansatz and the related gapless excitation spectrum implied by Eq. ( 15 



is in contradiction with the expected pseudogap physics at high temperatures [H] but 
this is a direct consequence of the choice of the ansatz. However, despite this deficit, 
the ansatz used here is in agreement with experiments that have been used for trying 
to observe the pseudogap phenomena O [54] . 

The other key approximation, namely neglecting the pairing part, is more critical 
in some cases than in others. Considering the values of the contact calculated in 
Section [6] as a benchmark, it appears that the full model in which also the imaginary 
part of the Hartree self-energy is included in the iteration works better in the case of 
the polaron than in the balanced system, whereas omitting the imaginary part works 
better in the balanced case than in the polaron case. This behavior is not at all 
surprising, considering that in the case of polarons, the many-body pair formation is 
strongly suppressed by the huge mismatch in the Fermi momenta. Hence the pairing 
part can be safely ignored and the Hartree self-energy contains all the relevant physics, 
including lifetime effects. Notice that the same argument holds also at sufficiently high 
temperatures when the pair formation is suppressed by thermal excitations, even in 
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the balanced gas. In the case of a balanced gas at low temperatures, the pairing part 
of the self-energy will also yield an imaginary part that will at least partly cancel the 
imaginary parts of the Hartree energy shift. In particular, close to the Fermi surface 
the cancellation is expected to be almost complete .34: . Thus neglecting the imaginary 
part of the Hartree self-energy in the case of a balanced gas is probably well justified. 

8. Conclusions 

In this manuscript I have calculated the Hartree self-energy of strongly interacting 
fcrmionic atoms. Short range of the atom-atom interaction potential and the diluteness 
of the gas allows an efficient description of the interactions through the momentum 
dependent scattering amplitude. I have calculated the Hartree energy shifts, effective 
mass, and the value of the contact for both balanced and spin-imbalanced gases. For 
the polaron system these values are in an excellent agreement with the results from 
other theories, including Monte Carlo results, and with experimental results, although 
it is not clear what value of the effective mass should be used as a reference. For the 
balanced system the agreement with more exact methods is also good when using 
Green's function ansatzes that satisfy the Bethe-Peierls condition. In the future, it 
will be interesting to consider the connection of the present theory with a BCS-type 
pairing theory by including also the pair part of the self-energy S. Furthermore, the 
generalization of the theory to inhomogeneous systems will be useful, allowing the 
treatment of the Hartree energy shift properly in Bogoliubov-deGennes-type spatially 
dependent mean-field theories (55] [56] . 
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Appendix A. Derivation of the many-body scattering T-matrix 

Here I will review the standard derivation of the many-body scattering T-matrix and 
how it can be expressed through the two-particle scattering amplitude. The calculation 
follows closely the derivation in Ref. [33] but here the discussion is kept more general 
as all the key approximations will be done in the main text. 

In the ladder approximation the many-body scattering T-matrix can be written 

as 



where Uq(P) is the interaction potential. Here Pi, P2 are the incoming four- momenta 
of the two particles and P^,Pa are the corresponding outgoing momenta (and while 
the above equation has been written without spin-indices, it is to be understood that 
Pi and P3 describe incoming and outgoing four-momenta of one spin state and P2 and 
P4 describe the other spin state) . This form for the many-body scattering T-matrix is 
not particularly useful since the potential Uq(P) is seldom known. However, for short 
range interaction potentials the T-matrix obtains a simple form when expressed using 
the scattering amplitude. 




x G(P 1 - Q)G(P 2 + Q)r(Pj - Q, P 2 + Q; P 3 , P 4 ), 



(A.l) 
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Appendix A.l. Scattering amplitude 

I will first review the standard scattering calculation of two particles of mass m and 
(relative) incoming momentum k scattering from a spherically symmetric potential 
Uo(r). I will restrict the analysis to the s-wave scattering, which is the relevant limit 
for low energy scatterings. 

The scattering of a particle, with incoming momentum k, from a spherically 
symmetric potential E/o( x ) is described by the Schrodinger equation 

^— (V 2 + k 2 ) </> k (x) = - (V 2 + k 2 ) V k (x) = C/ (x)^ k (x), (A.2) 

where x is the distance between the particles and the reduced mass m rc( j = m/2 is 
used in order to be able to describe the scattering of two particles with equal masses 
m. This equation can be rewritten as an integral equation 

V> k (x) = e* kx - J dy ?7 (y)Vk(y). (A.3) 

This can be expressed in the momentum space by taking the Fourier transform 

MP) = (27r) 3 <5(p - k) - h2 I . / ^ U (q)MP - q). (A.4) 

Defining the scattering amplitude 

- — /(s, k) = f A C/ (q)^k(s - q), (A.5) 

the scattering wavefunction becomes 

^(p) = (2nfS(p k) + - " ; lP ' J , , (A.6) 

2e, p - 2e fc - 277 

where e p = ^p 2 ■ The key quantity here is the scattering amplitude /(k', k) and the 
goal is to formulate the many-body scattering T-matrix using the scattering amplitude. 
However, the scattering T-matrix must yield the above scattering wavefunction for 
elastic scatterings and in the absence of many-body contributions. 

For future purposes, it is worthwhile to write down the completeness relation, 
which expresses the completeness of the scattering states in the absence of bound 
states, namely 

/ ^(P)^(P')* = ( 2 ^) 3(5 (P - Pi ( A - 7 ) 

The completeness relation, together with Eq. (A.5) yields 

I dk = _ 

J (ziry m 

Combined with Eq. (A.6 1 this yields for a hermitian potential Uq 

/(P,P')-/(P',P)* = - / 7^3 — /(P,k)/(p',k)* L - „ ^ r; 

J (27r) J m \2e k -2e p + ni 2e k - 2e p < - iij 

which can be considered as a generalization of the optical theorem for the inelastic 
scattering amplitude. 



,(A.9) 
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For short range potentials and elastic scatterings (k = s), the scattering amplitude 
is well approximated by 



/( S ,k)w— — — i -, forfc = S (A.10) 

1 + ika + ±r c fiak z 

where r e s is called the effective range and it is of the same order of magnitude 
as the interaction potential range. For a dilute atomic gas in the vicinity of a 
broad Feshbach resonance the factor fcr e g -C 1 and the effective range correction 
can be neglected. Indeed, for a zero-range delta function pseudopotcntial interaction 
Uo(r) — 4 ^ a ^( r )~3<r ( r ') tne scattering amplitude becomes 

/( 8j k)«— = °—, forfc = s. (A.ll) 
1 + ika 

However, while the scattering calculation is usually restricted to determining the 
elastic scattering amplitude, the inelastic scattering amplitude will be needed in 
the main text. I am not aware of a result for the inelastic scattering amplitude 
for the above delta function pseudopotential interaction, but I assume that the 
inelastic scattering amplitude depends only on the magnitude of the incoming relative 
momentum, i.e. 

/(s,k) = /(k)« TT ^-, Bar alls, k. (A.12) 
Such assumption for the inelastic scattering amplitude is consistent at unitarity a — > oo 



in the sense that it satisfies Eq. (A.9), but it is not necessarily true for any a and 



certainly not true for arbitrary potentials. 
Appendix A. 2. Many-body scattering T-matrix 



The T-matrix in Eq. ( A.l I can be expressed using a two-particle wavefunction Q as 
r(P 1) P 2 ;P 3 ,P 4 ) = J ^r i U (Q)Q(P 1 -Q,P i + Q;P 3 ,P 4 ), (A.13) 



which satisfies the ladder approximation for T in Eq. (A.l I if Q obeys 
Q(Pi,P 2 ;P 3 ,P4) = (2tt) 4 S(P 1 -P 3 )+iG(P 1 )G(P 2 ) 



J T^iWJfi^i-Q.i^ + Q;^,^), (A.14) 



which can be easily verified by substitution. The two-particle wavefunction is used 
here only in an auxiliary role and will soon be replaced by the pair susceptibility x- 

The two-body scattering processes are often easiest to describe using the center- 
of-mass coordinates. The total center-of-mass four-momentum satisfies 

2P = P 1 +P 2 =P 3 + P 4 , (A. 15) 

where the latter equality follows from the conservation of the total four-momentum 
in a homogeneous system (homogeneous also in time, hence the conservation of 
the frequency). The relative four-vectors are SP — |(Pi — P2) =: S and SP' = 
l(P3 — Pi) =■ K. The two-particle wavefunction expressed using these variables is 
now 

Q(P + S,P-S;P + K,P-K) = (2ir) 4 6(S - K) + iG(P + S) 

xG(P-S) J -0^U o (Q)Q(P + S-Q,P-S + Q;P + K,P-K) (A.16) 
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Assuming that the interaction potential Uq(Q) does not depend on the frequency 
(time-independent interaction potential), the frequency integral can be calculated in 
the right hand side, yielding 

Q(P + S,P-S;P + K,P-K) = (2tt) 4 5(S-K) 

+ iG{P + S)G(P -S) j C/o(q)x(5 - q, K; P). (A.17) 

where the pair susceptibility x is defined as 

X (S-c l ,K;P) = J ^ Q(P + S - Q,P - S + Q;P + K,P - K) 

= X (s-q,K;P). (A.18) 

The last equality follows from the observation that the pair susceptibility x does not 
depend on the frequency of the relative wave vector S. This is obvious if one does 
the change of integration variable Sq — qo — x. The equation does not depend on k 
either, as will become apparent soon. 

Substituting Eq. (A.17) into x yields 

X (s, K- P) = (2nf5(s - k) + i J ^ G(P + s - q )G(P -s + q ) 

X J ^(qMs-q^P). (A.19) 

This equation has no dependence on it is the defining equation for the function 
X, and the function definition remains exactly the same no matter what value is used 
for k . Choosing k — yields 



X (s, k; P) = (27r) 3 5(s - k) + J(s; P) J ^ U (d) X (s ~ q, k; P). (A.20) 
Not wanting to do key approximations in the Appendix, I have denoted here 

7(s; p) = l J^-G(P + s- qo )G(P - s + q ). (A.21) 
The scattering T-matrix can now be written in terms of x as 

r(p 1 ,p 2 ;P3,n) = J (Ss^oCq) J ~Q(Px-Q,P2 + Q;P 3 ,P4) 

= J Wf Ut){q) J ^-Q( p + S-Q,P-S + Q;P + K,P-K) 
= J ^o(q)x(s- q,k; P) 

= r(s,k;P). (A.22) 

Notice that the scattering T-matrix does not depend on the relative frequencies So 
and fcp because the pair susceptibility has no such dependence. 



Appendix A. 3. Two-body scattering T-matrix 

In order to obtain the connection to the scattering amplitude, one needs to first take 
the explicit two-particle limit of the many-body scattering T-matrix. In the case of the 
elastic scattering, this will yield the scattering amplitude. The two-body scattering 
T-matrix Tq is defined as the limit of T in which all the Green's functions G(K) are 
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replaced by the non-interacting Green's functions Gq(K) = l/(ko — £&) with the free 
particle dispersion e k = h 2 k 2 /2m. The two-body scattering T-matrix is thus defined 

as 



where 

Xo(s,k;P) = (27r) 3 5(s-k) 



r (s,k;P) = J -03[7 o ( q ) Xo ( s -q,k;P), 

dq 



(2tt) 3 
1 



(27T) 



t / o(q)xo(s - q, k; P), 



E P -2e 

and E P = 2p - Multiplying xo by (E P - 2e s ) gives 

(P P - 2e s ) Xo (s, k; P) - J £/ (q)xo(s - q, k; P) = 

(P P - 2e s )(27r) 3 (5(s - k) = (P P - 2e k )(2nf5(s - k 



(A.23) 
(A.24) 



(A.25) 

In order to obtain the connection to the scattering amplitude, the \o can be expressed 
as 

dk Vk(s)^(k)* 



Xo(s,k;P) = (Pp-2e fe ) 



(2tt) 3 



Ep - 2e- k 



(A.26) 



This can be confirmed by substituting this into Eq. (A.25) and remembering two 

i i 

1 



important equations: the Schrodinger equation in the momentum space (Eq. (A.4)) 

dq 



MP) = (27T) d 5( P -k)- 



2c, 



2e k -ir)J (2tt) 3 



C/o(q)^k(p-q),(A.27) 



and the completeness relation (Eq. (A.7)) 
dk 



(2tt) 3 



Mp)Mp'T = (27r) 3 <J(p-p'). 



(A.28) 



Substituting Eq. ( |A.26[ ) into Eq. ( |A.25| > yields 

^ ^(s)^(k)* 



(27r J )5(s-k) = (E P -2e s ) 



f/q 



(2tt) 3 
dk 



Uo(q) 



(2tt) 3 E P -2e~ k 

dk 7A £ (s-q)^ £ (k)* 



(2tt) 3 



E P - 2er 



(27T) 



(Pp-2e s )-^ S ^ 



E f 



2e E 



dq 



(27T) 



c / o(q)V'k(s- 



y> k (kr 

Pp - 2er 



dk 

(>]3 

(2e s - 2e s )^ k (p 



(Pp- 2£s )^ (S ^ (k) 



Pp - 2er 



> - 

kl °" kl "' [(Pp - 2e s ) - (2e fe - 2e s 



Pp - 2e s 

dk */> £ (s%(k) 
(2tt)3 Pp-2e^. 

dk 



7^k( s )^( k )* 
(27r) 3 <S(s-k). 



(A.29) 



Hartree shift in unitary Fermi gases 



28 



Eq. (A.26) provides now the required connection between the T-matrix 
formulation and the scattering calculation in |Appcndix A| There the solution of the 
Schrodinger equation Eq. ( A. 4 ) was connected to the scattering amplitude by Eq. ( A. 6 1 



tfk(p) = (2tt) 3 ,5(p - k) 



/(P,k) 

2e k - 2e p + in ' 



yielding now 



Xo(s,k;P) = -0 k (s)- 



Airh 2 



elk 



(27T) 



3^k( S ) 



1 



E P - 2e r 



2ei 



2e fc - in 



(A.30) 



/(k,k)*.(A.31) 



Similarly the two-body scattering T-matrix can now be expressed using the 
scattering amplitude as 



r (s,k;P) = J ^L_C/ ( q )V, u (s-q) 



Airh 2 
m 

Airh 2 
m 



(27T) 

/(s,k) 



C/o(q) 



dk 

(2^) 

2 



dk 



1 



Pf 



2er 



2<S 



/(s,k) 



26/c - i?7 
1 



P F 



2er 2e, - 2e k 



IT) 



(A.32) 
/(k,k)* 
/(k,k)*. 



A familiar form can be obtained when considering scatterings on-the-energy shell 
(scatterings that conserve energy). These correspond to the limit Ep = 2e k + irj, for 
which one obtains 

r (s,k;P) = /(s,k). (A.33) 



Appendix A. 4- Many-body scattering T-matrix revisited 

Here I will show the explicit connection between the many-body scattering T-matrix T 
and the two-body scattering T-matrix To, and hence also between T and the scattering 
amplitude /. This is of course important since the self-energy E is defined through T. 

The full many-body scattering T-matrix T was expressed using the pair 
susceptibility \ through Eq. (A.22). The \ can be rewritten as 

1 f dq 



X (s,k;P) 



Er 



2e< 



f/o(qMs-q,k;P) 



(2tt) 3 
= (27r) 3 <5(s - k) 



(A.34) 



/(s;P) 



1 



Ep — 2e s 



r(s,k;P) 



(2^) 3 «y(s-k) + x(s;P)r(s,k;P), 



(A.35) 

where the notation has been simplified by defining the regularized pair susceptibility 

1 



x(s;P)=7(s;P) 



Ep — 2e s 



(A.36) 



Recalling the definition of xo (Eq. (A. 24)), the operator part on the left hand side of 
Eq. (A.35 1 is equal to the inverse of xo, yielding the formal operator equation 

Xo k; P) X (k, k; P) = (2^) 3 <5(s - k) + x(s, k; P)T(s, k; P) (A.37) 

This can be solved using the Green's theorem and it gives 

x(s,k;P) = X o(s,k;P)+ J Xo (s, q; P) X (q; P)r(q, k; P)(A.38) 
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Multiplying by Uq(s — q) and integrating over q yields 

r(s,k;P) = r (s,k;P)+ J ^ro(s,q;P)x(q;P)r(q,k;P)(A.39) 

This equation provides the mapping between T and r . 

Assuming that the two-body scattering T-matrix ro(s,k;P) does not depend on 
cither of the two momenta s or k (but it may depend on the center-of-mass four- 
momentum P) allows writing the many-body scattering T-matrix in the standard 
form 

r(P) = i-r (P)/^x(q;P) ■ ( A ' 40) 
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